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Abstract 

The mixing of S — D partial waves for the heavy meson pairs in the decays T(5S) — > 
[B*B + c.c.]7r and T(5S) — > B*B*ir is considered. It is argued that this mixing, en- 
hanced by the heavy meson mass, is calculable as dominated by a rescattering through 
pion exchange if the production of the heavy mesons is dominated by iS-wave 'molecu- 
lar' resonances Z&(10610) and Z(,(10650). The effect of the mixing grows with the cm. 
energy E in each channel over the threshold, and may reach 10 -j- 20% in the rate at 
the upper end of applicability of the discussed approach, E ~ 15MeV. It is also argued 
that the mixing is likely to reach maximum at energy approximately equal to the mass 
gap between the thresholds. 



The strong interaction between hadrons containing heavy and light quarks results in a pe- 
culiar dynamics near a threshold for pairs of such hadrons. One of the recent most spectacular 
observations of this phenomenon is the discovery by Belle [1] of the 'twin' isotriplet reso- 
nances Zft(10610) and ^(10650) with masses within few MeV from the respective B*B and 
B*B* thresholds. The resonances were initially observed in the decays T(5S) — > T(nS)irir 
with n — 1,2,3 and T(5S) — > ht,(mP)7nr with m — 1 and 2 as conspicuous peaks in the 
invariant mass spectra for the system of bottomonium plus pion. These two peaks are nat- 
urally explained [2] as 'molecular' states with quantum numbers I G (J P ) = 1 + (1 + ) made 
out of S-wave pairs of mesons: Z b (10650) ~ B*B*, Z 6 (10610) ~ (B*B - BB*). This 
picture agrees well with the subsequent measurement [3] of the invariant mass distribution 
for the B*B + BB* and B*B* meson pairs in the decays T(5S) ->• [B*B + c.c.]^ and 
T(5S) — > [B* B*] ±r K T . The distribution in each channel displays a strong peak near the 
respective threshold, and in fact the data, although still with a large uncertainty, suggest 
that the entire distribution is dominated by the corresponding Z^ resonance. 

Naturally, more precise data, if available in the future, would allow a better analysis of 
the internal dynamics of the resonances. Such analysis however would require a better 
understanding of finer effects, including those in the spectra of the heavy meson pairs near 
the threshold resonances. The goal of this paper is to point out that the spectra in the 
vicinity of the resonances can be noticeably affected by an enhanced mixing of partial waves, 
specifically the S — D mixing, due to rescattering of the heavy mesons. In particular, in the 
(likely) case that the Zb resonances are indeed very strongly dominated by the .S-wave heavy 
meson pairs, a measurable .D-wave should arise due to rescattering at energies starting from 
the excitation energy above ~ 10 MeV over the threshold^. The onset of the D-wave in 
the spectrum is dominated by the pion exchange between the heavy mesons and is, to an 
extent, calculable. At higher excitation energy, however, the presence of the D-wave can 
only be estimated parametrically. The presence of the D wave can be measured in future 
experiments by studying the angular distribution of the heavy meson pairs and it is likely 
to be important for a better determination of the parameters of the Zb resonances. 

The enhancement of the effects of the interaction of heavy hadrons through the light 
degrees of freedom is well known. Indeed, in the limit of large mass M of the heavy quark 



1 The effect of the S — D mixing due to rescattering in the formation of the Zb resonances as 'molecular' 
bound states was recently considered in Refs. [4j [5] with the conclusion that this effect is small in the wave 
function of the meson pair in the bound state. 
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the latter interaction, described by a potential V, does not depend on M. Thus at a given 
momentum scale p any effect of the interaction enters through the product MV and thus 
gets bigger at large M. In particular, if V describes the spin-dependent interaction resulting 
in an S — D mixing, the effect of the mixing should generally contain the enhancement factor 
M/Aqcd in the amplitude. In a situation where the strength and the range of the potential 
are determined by Aqcd, an S-wave state rescatters into a D-wave state with the ratio of 
the amplitudes generally estimated as 

Ad Mp 2 

4- a A3 ' W 

A S A QCD 

where p is the cm. momentum of the heavy hadrons, and the estimate applies as long as 
p 2 <C Aq CD . This estimate, however, should be further modified for the effect of the pion 
exchange, since the pion mass fi can be considered as small in the scale of Aqcd- The 
contribution of the pion exchange to the wave mixing can be estimated as 



A D 2 M P ,2^2 ,,v, 

T K 9 t o at p < /i , (2) 
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and 
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— oc g - at A QCD > p > fi , (3) 
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where g is the constant for pion interaction with heavy hadrons. Clearly, these estimates 
imply that the pion exchange dominates the wave mixing at the energies corresponding to 
p 2 <C a qcd an d becomes comparable with the effect of other contributions to the spin- 
dependent interaction at the upper end of the applicability of both estimates (CQ) and (j3J) 
where the pion exchange can no longer be separated from those other short distance con- 
tributions. Neither the effect of the pion exchange at such momenta is calculable due to 
unknown form factors in both the pion interaction with heavy hadrons and in the short- 
distance behavior of As- For this reason the specific calculations in this paper are limited to 
the energies above the threshold corresponding to p 2 a qcd- Although restricted to this 
low energy range the calculation to be presented indicates that the arising from the rescat- 
tering D wave can be well measurable just above the peaks, corresponding to more than 
10% in the rate at p w 300 MeV, i.e. at the cm. energy E w 17MeV. Due to a moderate 
magnitude of the mixing one can neglect the back reaction, i.e. the feedback from the D 
wave to the S-wave, which approximation is assumed in the rest of this paper. 
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For the specific calculation in the discussed processes T(5S) — > [B*B + c.c.]7r and 
T(5S) — > B*B*7i one can write the pion interaction with the B and B* mesons in the 
form (see e.g. in Ref.[6]) 

+ ie ljk (Vjr a V k )} dm a (4) 

with Vi and P standing for the wave functions (in the nonrelativistic normalization) of the 
vector (B*) and the pseudoscalar (B) meson isotopic doublets and r a being the isospin Pauli 
matrices. The constant /„- ~ 132 MeV is used in Eq.(j4|) for normalization. The dimensionless 
pion coupling g can then be evaluated by using the heavy quark symmetry and the known [7] 
rate of the D* + — > Dix decay: g 2 « 0.18. 

The mixing of the partial waves due to the pion arising from the 'diagonal' rescattering 
in the channels B*B and B*B* is shown in the Figure 1. Only these 'diagonal' processes 
receive an unsuppressed contribution from the domain of the loop momentum q such that 
q 2 <C ^qcd provided that the overall cm. momentum p of the mesons satisfies the same 
condition. In particular, this is not the case for the rescattering between the two channels: 
B*B -H- B*B*, where due to the mass difference AM between the B* and B mesons the 
minimal momentum transfer through the pion propagator is of order q 2 ~ M AM, which 
corresponds to q 2 ~ ^qcd both numerically and parametrically. Thus the non-diagonal 
rescattering produces a mixing effect only of the order described by Eq.flTJ. 



H int = J- { (V^T a P) + h.C. 




Figure 1: The rescattering through the pion exchange resulting in an S — D mixing in the 
channels B*B and B*B*.) 



Using the notation a and b for the polarization amplitudes of the outgoing B* and B* 
mesons, one can write the general form of the T(5S) decay amplitudes: 



A[T(5S) ->• [B*B + c.c.]tt]/E^ 



V2 



As Sij + —j= A D (3 riiHj - 5ij] 
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A[T{5S) -> B*B*%}/E n 



V2 



A's $ij + A' D (3 riiiij — 5 i j / 



(5) 



where n = p/p is the unit vector in the direction of cm. momentum of one of the mesons, 
and p = \p\. The factor 1/E n in l.h.s. removes from the amplitude the chiral dependence on 
the energy E n of the emitted pion, which dependence is not the subject of the discussion in 
this paper. The amplitudes As and Ap are functions of p 2 and their relative normalization 
is chosen in such a way that the probability is proportional to the sum of their squares: 
\As\ 2 + |Ad| 2 - Under these notational conventions the amplitude of the D wave generated 
by the graphs of Fig. 1 can be calculated by means of nonrelativistic perturbation theory as 
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where in the latter transition an averaging over the relative angle between the vectors q 
and p is performed, so that only an integration over the scalar variable q 2 remains with the 
weight function F(p, q) given by 



F(p, q) 



5p 2 - 3q 2 - 3/i 2 4/i 2 p 2 + 3(p 2 - q 
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(7) 



In Eq.flE]) the notation A^ signifies that the relation can be used in either of the two heavy 
meson channels, and M stands for two times the reduced mass of the two mesons. For the 
purpose of the present calculation the difference in this reduced mass in the two channels 
can be neglected and a common value M w 5.3 GeV used. 

The absorptive part of the ratio of the D-wave and S-w&ve amplitudes, corresponding to 
the on-shell cut of the graphs of Fig. 1 can now be found unambiguously by setting q = p in 
the integrand in Eq.flS]): 



r = Im 
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Clearly, this expression exhibits the behavior described by Eqs.()2]) and ([3]) with one extra 
power of p/Aqcd arising from the phase space in the intermediate state. The plot of the 
ratio r given by Eq.fjH]) at g 2 = 0.18 is shown in Fig. 2 and illustrates the magnitude of the 
discussed effect of the S — D mixing. 
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Figure 2: The ratio r described by Eq.® as a function of the excitation energy above the 
threshold E = p 2 /M. 



One should notice that the expression (jSJ) does not describe all of the imaginary part 
of the ratio of the amplitudes and thus cannot be used to establish a lower bound on the 
D-w&ve generated through rescattering. The reason is that the amplitude A s \q 2 ) itself is 
complex with a phase generally depending on q 2 . In other words the graphs of Fig. 1 also 
generally have nontrivial on-shell cuts across the bubble depicting the 5-wave amplitude. In 
order to evaluate the full generated .D-wave amplitude a knowledge of the energy behavior 
of the S"-wave amplitude is needed. The data p] indicate that the Z& resonances dominate 
the spectrum of the heavy meson pairs in the discussed decays of T(5S). Thus it appears 
reasonable to use in Eq.(j5]) a pure resonance expression 

/ 2\ c const 

s{P ' = E-E + iY/2 = p 2 -pl + iMT/2 ' ( } 

with Eq being the energy position of the resonance relative to the threshold and Pq = 
ME . According to the current data the value of E is 2.7 ± 2.0 MeV for Z 6 (10610) 
and 1.8 ± 1.5 MeV for Z&(10650), while the corresponding values of T are 18.4 ± 2.4 MeV 
and 11.5 ± 2.2 MeV. Clearly, using the expression ((9]) for As makes the integral in Eq.(j6]) 
convergent and determined by a scale, which is a combination of p 2 , p^, MT/2 and fi 2 , all 
of which are assumed to be small in the scale of Aqcd and which justifies using the pion 
exchange amplitude as in Eq.([6]) without introducing a form factor. The integral can then be 
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readily calculated, and the plots in Fig. 3 illustrate the magnitude of the discussed generated 
D-wave. For this illustration the value E = is used and the effect depends very weakly on 
this particular value as long as £0 is smaller than T/2. Also for this illustration the value 
T = 14MeV is used, which is close to the data for both resonances. One can see that 
at the upper end of applicability of the discussed approach, i.e. at E ~ 15 -j- 17MeV the 
relative contribution of the generated D wave can amount to 10 -j- 20% in the event rate. 



D/S 




E 



Figure 3: The absolute value (solid) and the real (dashed) and imaginary (dotted) parts 
of the ratio Ap/Ag for the D-wave amplitude generated from the As given by Eq.(|9]) with 
E = and T = 14 MeV. 

It can be noticed that the simplest Breit-Wigner formula for the resonance amplitude is 
used for As in Eq.([9]). In particular the width parameter T is assumed to be constant and 
does not include the energy- dependent absorption in the heavy meson channel above the 
threshold. This approximation appears to be sufficient with current data, but may need to 
be modified in a thorough analysis of possible more detailed data in the future. 

As one can see from the plots of Fig. 3, the S — D mixing is quite small at low energies 
E ~ 1 MeV, which are relevant for the internal dynamics of the resonances, so that 
the mixing is likely only a minor effect in this internal dynamics, in agreement with the 
conclusions of Refs. jU |5]. The ratio Ap/As however grows with the excitation energy and 
can possibly be studied in the decays from T(55). At present there seem to be no theoretical 
means of analyzing the behavior of the S — D mixing at larger excitation energies E beyond 
the region of applicability of the presented approach, i.e. at E being not small as compared 
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to Aq CD /M. It is however interesting to notice that if higher values of E were accessible in 
the T(5S) decay, i.e. in an artificial theoretical limit of a large mass of T(5S), the S — D 
mixing should vanish at high E. Indeed, the strong interaction rescattering arises from the 
interaction of the light quarks in the heavy mesons, and the spin of the heavy b quarks is 
conserved, and the whole process can be considered with spinless b quarks. With spinless b 
quarks the T(55) would be an J p = + state, and after emission of an S- wave pion would 
produce heavy meson pair with quantum numbers J p = 0~, for which obviously no S — D 
mixing is possible. In terms of the mesons, the channels B*B and B*B* would degenerate 
into one channel as well as the two Z\> resonances would coalesce into one. The vanishing 
of the S — D mixing then occurs through a cancellation between the 'diagonal' rescattering, 
e.g. B*B* — > B*B* and the 'off-diagonal' B*B — > B*B*. For the actual heavy mesons, 
whose thresholds are split by AM, in the region of dominance of the pion exchange only 
the 'diagonal' scattering should be retained, as discussed previously. However at large E, 
specifically at E 3> AM the difference in the thresholds becomes unimportant and the heavy 
quark spin symmetry behavior should set in. Therefore one can reasonably expect that the 
D/S ratio has a maximum at E ~ AM and decreases at higher energy. It would be quite 
interesting if this behavior can be studied in the decays from the actual T(5S). 

It can be also mentioned that a similar enhanced mixing of partial waves, namely the 
P — F mixing, can be expected in the production of B*B* pairs in e + e~ annihilation at the 
cm. energy just above the threshold. A detailed discussion of the effects of this mixing will 
be reported separately. 
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